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Abstract. We consider an orthonormal basis of eigenfunctions of the Dirichlet Laplacian 
for a rational polygon. The modulus squared of the eigenfunctions defines a sequence of 
probability measures. We prove that this sequence contains a density-one subsequence that 
converges to Lebesgue measure. 



1. Introduction 



One of the key challenges in quantum chaos is to understand how quantum eigenstates 
distribute in the semiclassical limit. In the present note we will discuss this problem for 
eigenfunctions of the Dirichlet Laplacian for rational polygons in the plane, with particular 
attention to the work of Bogomolny and Schmit on "superscars" [4J. The conclusion of the 
study presented here is that any such scarring may, for almost all eigenstates, only appear in 
momentum, but not in configuration space. 
1^ ' Let Z? C be a bounded connected domain with piecewise smooth boundary. The Dirich- 

. let Laplacian A/) is defined as the standard Laplacian acting on functions in C^{D) that vanish 

at the boundary of D. The eigenvalues of the positive definite operator — A^i (which plays 
the role of the quantum Hamiltonian) will be denoted hy < Ei < E2 < < ■ ■ ■ ^ 00, and 
^ , the corresponding eigenfunctions by ipi, -02; V'S) etc. A striking result on the asymptotic dis- 

I tribution of eigenfunctions is the Schnirelman-Zelditch-Colin de Verdiere quantum ergodicity 

' theorem, which in the present setting is due to Zelditch and Zworski }15] : If the billiard flow 

ly-^ i is ergodic in S*D (the unit cotangent bundle of D), then almost all eigenfunctions become 

' uniformly distributed in S*D. 

Hassell [5j has shown that for certain domains D with ergodic flows there are subsequences 
of eigenfunctions which fail to become uniformly distributed. Thus the restriction to subse- 
quences in the above quantum ergodicity theorem is in general necessary. We are not aware 
of examples of domains D for which the full sequence of eigenfunctions becomes uniformly 
distributed. 

^ ■ In the present note we point out that uniform distribution of the eigenfunctions in configu- 

ration space (rather than the full phase space) holds for domains given by rational polygons. 
A rational polygon is a simple plane polygon (which means that its interior is connected and 
simply connected), so that all the vertex angles are rational multiples of it. Billiards in rational 
polygons give rise to dynamics that is neither integrable nor ergodic on S*D (except of course 
for the few integrable cases — rectangles, equilateral triangle, and right triangles with an acute 
vertex angle of either 7r/3 or vr/4). The absence of ergodicity is due to the fact that, for any 
initial direction, the motion in S*D is to a higher genus flat surface. By exploiting the ergodic 
properties of directional flows on such surfaces, we can however show that quantum ergodicity 
still holds in configuration space: 
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Theorem 1. Assume D is a rational polygon. Let {ipn} orthonormal basis of eigenfunc- 
tions of the Dirichlet Laplacian on D. Then there is a density- one sequence rij S N such that 
for any subset A C D with boundary of measure zero, 

(1) lim / |V^„^,(x)pdx-^'"^(^) 



area(D) 



Density-one means here that 



(2) Urn 



The proof of Theorem [T] follows Zelditch and Zworski's approach in |15j until the last step, 
which is the only place where ergodicity of the billiard flow on S*D is used. Up to that point 
all arguments hold for any pseudo-differential operator of order zero. We then specialize to 
multiplication operators and appeal to the theorem of Kerckhoff, Masur and Smillie [8j who 
showed that for rational polygons, almost all directional flows (see ^ are uniquely ergodic. 
Details of the proof of Theorem [1] are provided in ^ 

The study of polygonal billiards in the context of quantum chaos goes back to Richens and 
Berry in [10], who understood that the "pseudo-integrable" nature of the classical dynamics 
has an important effect on the quantum spectrum. The energy level statistics appear to be 
intermediate between those believed valid for generic chaotic and integrable systems, in that 
the level spacing distribution shows level repulsion at small distances and Poisson tails at 
large distances [l0l[T2l[3]. The quantum wave functions were investigated empirically, both 
numerically [HITS] and experimentally using microwave cavities [9l[2], finding "strong scarring" 
related to families of periodic orbits. Bogomolny and Schmit ^ constructed long-lived states 
(quasimodes) associated to families of periodic orbits, which they called "superscars" , and 
suggested that a positive proportion of true eigenfunctions have large overlaps with such states 
at high energies. Our Theorem [1] shows that this phenomenon cannot occur in configuration 
space for subsequences of density larger than zero. 

A recent rigorous result that holds for all eigenfunctions is a bound by Hassell, Hillairet 
and Marzuola [6] which establishes that eigenfunctions cannot localize their mass away from 
the polygon's vertices. 



Remarks. 

(i) Theorem [1] seems new even for the classical integrable case of the square billiard. In 
that case a proof based on harmonic analysis and simple arithmetic considerations is 
also available. See [7| for information on the possible quantum limits of exceptional 
subsequence of eigenfunctions. 

(ii) We can relax the requirement in our theorem that the polygon D is simple, i.e., that 
its interior is simply connected. (We still require that D is connected.) Then the 
correct definition of a rational polygon is that the group generated by the linear parts 
of the reflections in the sides of the polygon is finite. This implies that all vertex angles 
are rational multiples of vr, and is equivalent to the definition in the simply-connected 
case, but not necessarily in the multiply connected case. 

(iii) Theorem [T] also holds for the Neumann Laplacian of D, and more generally for the 
Laplacian of an arbitrary translation surface. 

(iv) A result of the same nature was recently obtained in [HI, by completely different 
methods, for a point scatterer on the torus. For this system, the classical dynamics is 
integrable, but the quantum problem is not. 
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2. Background on billiards in rational polygons 

For billiards in a rational polygon D, the phase space is the unit cotangent bundle S*D, 
which is a direct product 



(3) 



S*D = DxS^ = {{x,u 



The normalized Liouville measure is defined as 



(4) 



area(D) 



dx d(j). 



The billiard fiow is defined on S*D via specular refiection for all trajectories not hitting the 
vertices of the polygon D; the refiection law for trajectories that hit a vertex can be defined 
arbitrarily. (The latter trajectories form a set of measure zero and are ignored in the following 
discussion.) The measure d/i is invariant under the billiard fiow. 

Let r be the group generated by the linear parts of the refiections in the sides of the polygon 
L), which is a finite group since our polygon is simple and all vertex angles are rational multiples 
of vr. For each direction 0, the set 



(5) 



De:=Dx |J {7^} 
76r 



is preserved by the fiow; call this restriction $g the directional flow. Kerckhoff, Masur and 
Smillie showed |8j that for almost all 9, the directional fiow <I>g is uniquely ergodic. 

Observables are (smooth) functions a{x,uj) on S*D. "Isotropic observables" are observables 
which depend only on the position variable x, that is are independent of the momentum 
(direction): 

(6) a(x, bj) = a^ix). 

The time average of an observable a G C^{S*D) is 



(7) 



1 f'^ 

a^ix,u}):=— / ao ^*{x,uj)dt. 

2T J_rp 



Lemma 2. Let a{x,uj) = ao(x) be an isotropic observable. Then the time averages satisfy 

2 



(8) 



lim 



|ar| dfi 



1 



area(L') Jd 



aQ[x')dx' 



Proof. As a consequence of the Kerckhoff-Masur-Smillie theorem, there is a set of full Lebesgue 
measure Q. <Z S^, such that for all (x, uj) G D xVt, the time averages a"^ converge to the average 
of a on D^: 



(9) 



lim a (x, oj) 

T^oo 



adfi^ = V , 

^ area(D) Jo 



a{x' , ^oj)dx' . 



By the dominated convergence theorem we find that 



(10) 



lim 



\aT\ dfj. 



S*D 



#r ^ area(Z)) ^ 



a{x' , ^uj)dx 



doj. 



We conclude by noting that for isotropic observables, (fTOj) reduces to 



□ 
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3. Proof of Theorem [T] 

We assume throughout this section that D is a rational polygon, and {ipn} is a fixed 
orthonormal basis of eigenfunctions of the Dirichlet Laplacian on D. 

By a standard density argument, Theorem [1] is a direct consequence of the following variant 
for smooth isotropic observables a{x,uj) = ao(x) with compact support in D° (the interior of 
D): 

Theorem 3. There is a density- one sequence uj G N such that for any oq G C'^{D°), 

(11) lim / ao{x)\i}nAx)\'^ dx = Kr^l ao{x)dx. 

J^'^Jd aiea{D) Jjj 

Let ^ be a pseudo-differential operator of order zero with principal symbol a. The quantum 
variance for A relative to the basis {ipn} is defined as 

(12) V{A,E) := \ |(^V'n,V'n> - a|', 

with a := Jg,^a{x,oj) dfi^x^oj) being the phase-space average of the observable a. We recall 
the local Weyl law in this context \15\ Lemma 4]: 

(13) jim ^p;^ ^5^^^^"' ^'^^ = 
As explained in [15] , to prove Theorem [3l it suffices to show 

Theorem 4. Let A be the multiplication operator Af{x) = ao(a;)/(x) with gq G C^(D°). 
Then 

(14) lim V{A,E) = 0. 

Proof. We follow the proof of quantum ergodicity for billiards given in pT^. (This elegant 
argument was first developed in [14j.) The key bound we require is 

(15) limsupF(^,^) < / -afdi^ 

E^oo Js*D 



for any T > 0. Eq. (|15p holds for any billiard with piecewise smooth boundary, and any 
pseudo-differential operator A whose Schwartz kernel is compactly supported in D° x D°. 
The proof of (jlSp follows directly from the estimates in [15] ; we will summarize the main steps 
below. Crucially, (jlSp does not require any assumptions on the ergodicity of the billiard flow. 
We now specialize to A whose principal symbol is an isotropic observable a{x,uj) = ao{x), 
and input the (unique) ergodicity of almost all the directional flows in the form of Lemma [2j 
Thus, for any e > 0, there exists T = > so that fg,jy — a^d/i < e and therefore 
lim.supE^^V{A, E) < e. This proves Theorem [H □ 



Outline of the proof of eq. p5|) : Fix T > 0. Let Us C S*D be the set of initial conditions 
that, within the time interval [— 2T, 2T], do not hit a ^-neighbourhood of the vertices of the 
polygon D. We decompose A = As + Rs, where ^4^ is chosen so that (i) its symbol has essential 
support in Us, and (ii) the measure of the essential support of the symbol of Rs tends to zero 
as (5 — )■ 0. The principal symbol of As is denoted by as- Since ^pn is an eigenfunction of Ad 
we have 

^^g^ K(^<5 - a<5)V'n,'0n>P = K(^<5 - a5)T^n,^n>P 

< {{As -as)T{As -as)T'>Pn,fpn), 

where {As)t '■= ^ J^rp ex.p{it\/—A[))As exp{—it\/—Ar))dt. Egorov's theorem [TJl Lemma 5] 
implies that the principal symbol of {As)t is aj. This, in conjunction with the local Weyl 
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law ([T3]) applied to {As — as)T{As — as)^, proves eq. p5|) for the truncated A^. The result is 
extended to the original A by showing that both sides of the inequality are obtained as the 
5 — )■ limit of the truncated version. The central ingredient in the required estimates is again 
the local Weyl law; see [15] for details. 
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